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Abstract 



We report a measurement of the longitudinal alignment Al of the recoil 
nucleus in the /x"+ 14 N(l+, 0) -> i/ M + 14 C(2 + , 7012) transition. The exper- 
iment was performed on the M9B beamline at the TRIUMF cyclotron via 
the measurement of the Doppler lineshape of the subsequent 14 C(2 + , 7012) — > 
14 C(0 + ,0) gamma-rays. We compare our result Al = 0.60 ± 0.11 to various 
model calculations, and discuss the sensitivity to the induced pseudoscalar 
coupling, second-forbidden effects and 2hu> wavefunction admixtures. 
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I. INTRODUCTION 



Muon capture by complex nuclei has been used in the study of both the dynamics of the 
weak interaction and the structure of the atomic nucleus. Generally, in exclusive capture 
between discrete states one attempts to separate the effects of dynamics and structure by 
suitable choice of physical observables and spin-parity sequences. Of course in practice these 
dual aspects of muon capture are not completely separable, and inevitably such studies must 
confront the entanglement of dynamics and structure. 

Of special interest in muon capture is the induced pseudoscalar coupling constant g p of 
the nucleon's weak axial current A^. Its value is firmly predicted by chiral symmetry [1-3] 
thus making its measurement an important test of low energy QCD. Moreover, in complex 
nuclei its medium modification - due to effects that range from core polarization and ex- 
change currents to partial restoration of chiral symmetry [5-7] - are of considerable interest. 
Unfortunately, the world data on g p from exclusive capture on complex nuclei is rather 
sparse, it comprising the study of recoil orientations in allowed Gamow- Teller transitions on 
12 C [8-11] and 28 Si [12-14], the study of hyperfine effects in allowed Gamow- Teller transi- 
tions on n B [15,16] and 23 Na [17,18], and the measurement of the 16 O(0 + ,0) ^ 16 N((T, 120) 
first-forbidden transition rate [19-21]. While the majority of these experiments are in agree- 
ment with the chiral prediction for the induced coupling, the result of g p /g a = l-O^ii [4] 
from the 28 Si experiments is in stark disagreement with the theoretical expectation. 

It is perhaps surprising - given that fifty years have passed since the birth of V- 
A theory - that few measurements of recoil orientations in muon capture have actually 
been performed. Partly motivated by the puzzling result for the recoil orientation in the 
28 Si(0 + ,0) ^ 28 A1(1 + , 2201) experiment, we therefore decided a further investigation of re- 
coil orientations in exclusive capture was worthwhile. The case we chose was 14 N(l + ,0) 
— > 14 C(2 + , 7012), it presenting an interesting example of the different sensitivities of the 
different observables to the weak dynamics and the nuclear structure. 

Herein we report the measurement of the recoil alignment in the transition 14 N(l + ,0) 
— > 14 C(2 + , 7012). In Sec. II we briefly outline the measurement technique and experimental 
setup. The determination of the recoil alignment from the Doppler broadened spectrum 
of the 7012 keV gamma-rays is described in detail in Sec. III. In Sec. IV we compare 
our experimental result to model calculations, and discuss the sensitivity to the induced 
pseudoscalar coupling g p , second-forbidden contributions, and A = 14 nuclear wavef unctions. 
Note that gamma-ray intensities from this /z _14 N experiment were published earlier in Stocki 
et al. [22]. 

II. DOPPLER METHOD AND EXPERIMENTAL SETUP 

The products of exclusive muon capture are a left-handed muon neutrino and an orien- 
tated recoil nucleus. The recoil orientation is a direct manifestation of the V-A character 
of the weak interaction, and for certain spin-parity sequences the induced coupling g p has a 
large effect on the recoil orientation. Moreover, in appropriate cases the recoil orientation 
can be experimentally determined from the directional correlation that is imparted on either 
the beta-rays or the gamma-rays that are subsequently emitted by the unstable recoil. 
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A method for measuring the 7-recoil directional correlation was originally proposed by 
Grenacs et al. [23]. The method takes advantage of the Doppler shift of the 7-ray energy 
for a decay in-flight. If the recoil is in motion as it decays its energy is shifted by 



AE E-E Q 



= (3 cos 9 



(1) 



E Q E Q 



where E Q is the 7-ray energy in the recoil reference frame, E is the 7-ray energy in the 
laboratory reference frame, f3 = v/c is the velocity of the recoil in the laboratory, and 9 is 
the angle between the 7-ray momentum vector and the recoil momentum vector. Conse- 
quently the lineshape of the Doppler broadened 7-ray is a reflection of the 7-recoil angular 
correlation. In suitable cases, i.e. where the 7-ray lifetime is short and the spin-parity se- 
quence is sensitive, this method thus permits the determination of the recoil orientation in 
the capture process. 

In this work we have studied the sequence 14 N(l+,0) -> 14 C(2+,7012) -> 14 C(0+,0). 
Because the lifetime of the 7012 keV state is rather short, 9.0 ± 1.4 fs [24], and the slowing- 
down time of the 14 C recoil in a liquid N 2 target is rather long, ~0.8 ps [25], the 7-ray 
spectrum is Doppler broadened. The resulting 7-recoil angular distribution is given by 



where P 2 (cos9) is the Legendre polynomial, Al is the recoil's longitudinal alignment and 
B 2 is the 7-decay correlation coefficient. 1 Note that the longitudinal alignment Al is a 
direct manifestation of the different populations of the magnetic sub-states about the recoil 
direction, and is entirely governed by the capture process. Conversely the 7-decay coefficient 
B 2 is a function of the 7-decay. For a 2 + — > + transition and a rank-2 orientation its value 

is B 2 = — \ ji [29]. For further details see Appendix A. 



A detailed discussion of the experimental setup was presented in Stocki et al. [22] - herein 
we review the features that are important for the determination of the recoil alignment Al 
in the 14 N(l + ,0) ^ 14 C(2 + , 7012) transition. The experiment was conducted on the M9B 
superconducting muon channel at the TRIUMF cyclotron. The negative muon beam had a 
momentum of 65 MeV/c and a flux of 2 x 10 5 s -1 , with an electron contamination of ~20% 
and a pion contamination of <0.2%. The incoming particles were detected in a three-element 
plastic scintillator beam telescope and stopped in a liquid nitrogen target. The target was 
20 x 55 cm in cross section and 21 cm in thickness with walls of polystyrene. Mu-metal 
was used to reduce the ambient magnetic field and thereby the muon spin precession. The 
outgoing 7-rays were detected at right-angles to the beam axis in a high-purity Ge detector 
with a Compton suppressor. The Ge detector had an efficiency of 44%, an in-beam energy 
resolution of 2.5 keV, and an in-beam time resolution of 6 ns, for the 1.33 MeV 7-ray line 
of a Co-60 source. The Compton suppressor comprised a six-element annulus of Nal(Tl) 
crystals. 



In Eqn. 2 we omit the negligible effects of rank-4 correlations and assume either unpolarized 
muonic atoms or perpendicular geometry (see Appendix A and Ref. [4] for further details). 



[26-29] 



W{9) = l + A L B 2 P 2 (cos9) 



(2) 
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A remark is warranted on the hyperfine effect in the /i _14 N atom. On formation the IS 
ground state of the /j,~ 14 N atom is a statistical mix of the two hyperfine states, i.e. one 
third F_ = 1/2 and two thirds F + = 3/2. Hyperfine transitions, from the higher-lying 
F + state to the lower-lying F_ state, would modify these proportions, and consequently 
the various observables in muon capture. Such transitions in //-atoms are mediated by 
the Ml emission of an Auger electron from the surrounding electronic orbits [30,31]. For 
yU _14 N, where the muonic atom hyperfine splitting is 7.4 eV and outer-most electron binding 
energy is 11.3 eV, one expects such transitions to be energetically forbidden [22]. We note 
the authors of Ref. [32] have suggested the possibility of a non-zero hyperfine transition 
rate via an unconventional transition mechanism in the /i~ 14 N atom. However, in absence 
of confirmation in recent studies by Stocki et al [22], we will hereafter assume that the 
capture originates from a statistical mixture of the hyperfine states, i. e. the measured recoil 
alignment is the statistical recoil alignment A s ^ at . 

Finally we note the proximity of the interesting 7012 keV level to the 14 C neutron sepa- 
ration energy. This allays the concern over cascade feeding to the 7012 keV level from muon 
capture to a higher-lying level. 

III. ALIGNMENT DETERMINATION 

The resulting spectrum in the energy region of the 7012 keV 7-ray is shown in Fig. 1. 
The 7-ray events in Fig. 1 were recorded within a ~ 5.0 /is time window of a muon arrival 
and in absence of any signal in the Compton suppressor. 2 

The first step in the analysis was the determination of the energy calibration and the 
instrumental resolution of the Ge detector. The energy calibration was obtained from eight 
well-known and clean 7-ray lines from thermal neutron capture on Al, Fe and CI with 
energies that ranged from 6.2 MeV to 7.8 MeV [33]. In fitting the peaks we found a simple 
Gaussian lineshape with energy-independent width to be sufficient. A linear relation between 
channel number and 7-ray energy was also adequate for the calibration in the energy range 
6.2-7.8 MeV. 

Next we performed the least-squares fit to the Doppler lineshape of the 7012 keV 7- 
ray. The fit function involved the convolution of the theoretical lineshape (Eqn. 2) with 
the instrumental resolution of the Ge detector and the slowing-down effects in the target 
material. To compute the slowing-down time r s d of the 14 C recoil ion in the liquid N 2 
target we used the computer program SRIM [25]. We note that the slowing-down time 
t S( i ~ 0.8 ps is much longer than the gamma-ray lifetime r = 9.0 ± 1.4 fs, and consequently 
in incorporating the slowing-down effects into the Doppler lineshape we assumed a simple 
linear form (3{t) = f3 Q (l — t/T sd ) for the recoil velocity time dependence. 

The energy spectrum in Fig. 1 also indicates the presence of a weak background line on 



2 Note that the first escape peak and the second escape peak of the 7012 keV gamma-ray were also 
identified, and therefore in principle were available for the determination of the Doppler lineshape. 
However, in practice such analysis was precluded by the poor statistics and the line backgrounds 
in the regions of these peaks. 
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the far lefthand side of the Doppler spectrum. We identified this peak as a 6978 keV gamma- 
ray originating from thermal neutron capture on chlorine-35. This identification was based 
on both its energy correspondence and reasonable consistency with observed intensities of 
other stronger 35 Cl(n, 7) lines. The 35 Cl(n, 7) background was believed to originate from 
polyvinyl chloride tape in the vicinity of the target. Consequently, both a Doppler peak and 
a Gaussian peak were needed to fit the region of the 7012 keV peak. 

The fit function involved a total of eleven adjustable parameters. They comprised the 
amplitudes (A D , A G ), centroids (X G , X D ) and instrumental widths (a) of the Doppler peak 
and the Gaussian peak, the initial velocity slowing-down time (r sd ) and longitudinal 
alignment (A s L tat ) of the 14 C recoil, and the amplitude and the slope of the continuum 
background. Fortunately, a number of parameters could be fixed in the fit, for example the 
instrumental widths and peak centroids from the energy calibration, the initial velocity from 
the reaction kinematics, and the recoil slowing-down time from the SRIM calculation. Also 
the amplitude of the 6978 keV 35 Cl(n, 7) line could be estimated from the amplitudes of the 
other 35 Cl(n, 7) lines and the available thermal neutron capture data [33,34]. 

Our benchmark fit, which clearly demonstrates the good agreement between the mea- 
sured spectrum and the theoretical function, is shown in Fig. 1. In this fit the instrumental 
width (a), initial velocity and slowing-down time ij sd ) were all fixed at their calculated 
values. The fit gave A s L tat = 0.60 ± 0.11 with a chi-squared x 2 p df = 0-7- 

To investigate the correlations between the recoil alignment and the other parameters 
we performed a series of fits, the results of which are summarized in Table I. The Table 
shows that the "best-fit" value of the recoil alignment is only weakly dependent on the 
parameters related to the energy calibration and the instrumental resolution (rows 1-4 of 
Table I). In addition it shows that large changes in the input value of the slowing-down 
time, which ranged from r sd = 0.4 ps to r sd — > 00, had very little effect on the "best-fit" 
value of the recoil alignment (rows 5-6 of Table I). Moreover, the choice of either a flat 
continuum background or a linear continuum background and the left-right margins of the 
fit region, had only minor effects on the "best-fit" value of the recoil alignment. 

By contrast the presence of the 6978 keV background line had a significant impact on the 
least-squares-fit to the recoil alignment. The importance of including the 6978 keV line is 
demonstrated by comparing the results of fits performed with the 6978 keV peak included, 
which yielded A L tat = 0.604 ± 0.106 (row 1), and the 6978 keV peak omitted, which yielded 
j^stat _ Q_42g -j- 0.078 (row 8). In addition, we found a significant difference in "best fit" 
values of A s L tat when comparing the cases of either (i) varying the 6978 keV line amplitude 
or (ii) fixing the 6978 keV line amplitude. In order to fix the amplitude of the 6978 keV 
35 Cl(n, 7) line we used the amplitudes of the nearby 6678 keV and 7414 keV 35 Cl(n, 7) lines 
and their relative intensities taken from Krusche et al [34]. The free-amplitude fit gave a 
recoil alignment Al = 0.604 ± 0.106 with a 6978 keV peak amplitude A G = 89 ± 30 (row 
1). By comparison the fixed-amplitude fit gave A L = 0.503 ± 0.083 with A G = 41 (row 9), 
i.e. showing a difference of ~1 a for the recoil alignment and a discrepancy of ~1.6 a for 
the peak amplitude. 

The 1.6 a discrepancy between the estimated amplitude and the fitted amplitude of the 
6978 keV line may reflect the statistics of the 6978 keV background line, the presence of an 
unidentified background line, or the uncertainties in interpolating 35 Cl(n, 7) intensities. We 
consider the result of the free-amplitude fit to be most trustworthy. 
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In conclusion, by fitting the 7012 keV lineshape we obtained the recoil alignment A s L tat = 
0.60 ± 0.11, which indicates the 14 C(2 + , 7012) recoil is highly orientated following muon 
capture. We found the alignment determination was relatively insensitive to the parameters 
related to the energy calibration, instrumental resolution, and slowing-down effects, but 
more sensitive to the presence of the 6978 keV 35 Cl(n, 7) background line. 



IV. INTERPRETATION 



A. Theoretical overview 



The mass- 14 system includes a well-known example of a highly suppressed allowed 
Gamow- Teller transition; the /5-decay that connects the 14 C(0 + , 1) ground state and 
14 N(l + ,0) ground state. In lieu of g.s.-to-g.s. Gamow- Teller strength the 14 N^ 14 C GT 
strength is distributed across a few low-lying 14 C excited states, with most notable the 2 + 
doublet at 7012 and 8318 keV [35]. Interestingly, the 2 + doublet wavefunctions involve both 
(Op)' 2 components and (0p)~ 4 (ls0d) 2 components [36-38], with configuration mixing being 
indispensable to understanding the partitioning of the GT strength within the 2 + doublet. 

Calculations of the /i"+ 14 N(l + , 0) -> i/„+ 14 C(2 + , 7012) capture rate have been performed 
by Bukhvostov et al. [43], Mukhopadhyay [44,45], Kissener et al. [46], Desgrolard et al. [47,48] 
and Auerbach and Brown [49]. It is clearly established that a simple (Op) -2 model for the 

+ 14 iV(l + , 0) — > f (U + 14 C(2+, 7012) transition over predicts the rate by roughly a factor of 
5-10. As alluded to above, this discrepancy is understood as a result of the importance of the 
(0p)~ 4 (ls0d) 2 components in the 14 C(2 + , 7012) wavefunction, such configurations decreasing 
the wavefunction overlap between the 14 N initial state and the 14 C final state. 

The longitudinal alignment of the 14 C recoil in the 14 N(1 + , 0)^ 14 C(2 + , 7012) transition 
has been studied theoretically by Bukhvostov and Popov [26,27] and Dmitrieva et al. [28]. 3 
Bukhvostov and Popov [26,27] were first to advocate the alignment as a possible probe of 
the proton's induced pseudoscalar coupling g p . Dmitrieva et al. [28] then extended this 
work to calculate both the recoil alignment and its hyperfine dependences using a (Op) 2 
model space and including the forbidden nuclear matrix elements. For g p /g a ~ 7 they 
obtained an alignment A s L tat ~ 0.60 when including the second-forbidden corrections and an 
alignment Af at ~ 0.23 when excluding the second-forbidden corrections, the authors tracing 
the striking contribution of forbidden terms to the operator M 2 2 • cr. The possible effects of 
the (0p)~ 4 (ls0d) 2 components of the 14 C(2 + , 7012) wavefunction on the 14 C recoil alignment 
were however not considered in Refs. [26-28]. 



3 Actually the authors give results for the 7-recoil angular correlation coefficient in the 14 N(1 + , 0) 
— ► 14 C(2 + ,7012) — ► 14 C(0 + ,0) sequence. However, as discussed in Appendix A, the 7-recoil corre- 
lation coefficient 02 and recoil alignment Al are related according to 02 = ~\f^A^L for 14 N(1 + , 0) 
-> 14 C(2+,7012) -> 14 C(0+,0). 
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B. Model calculations 



Herein we have performed shell model calculations of the capture rate and the recoil 
alignment for the 14 N(1 + , 0)^ 14 C(2 + , 7012) transition in both a (Op) -2 model space and a 
(0p)- 2 +(0p)- 4 (ls0d) 2 model space. For the (Op)" 2 calculation (denoted CKPOT) we used 
the Op-shell effective interaction of Cohen and Kurath [50]. For the (0p) _2 +(0p)~ 4 (ls0d) 2 
calculation we used the two-body matrix elements of: Cohen and Kurath [50] for the Op 
shell interaction, Chung and Wildenthal [51] for the IsOd shell interaction, and Millener and 
Kurath [52] for the cross-shell interactions (for further details see Warburton and Millener 
[53]). Moreover, we performed three versions of the (0p)~ 2 +(0p)~ 4 (ls0d) 2 calculation. In 
version one (denoted MK3CW1) the IsOd single particle energies were fixed as given in Ref. 
[53]. In version two (denoted MK3CW2) the ls-Od single particle energies were shifted by 
—3.0 MeV (from Ref. [53]) in order to reproduce the ~50% admixture of the (0p)~ 4 (ls0d) 2 
configurations in the 14 C(2 + ,7012) wavefunction as determined in Refs. [36-38]. In version 
three (denoted MK3CW3) the ls-Od single particle energies were shifted by —4.9 MeV 
(from Ref. [53]) in order to reproduce the 14 N(1 + , 0)^ 14 C(2 + , 7012) capture rate of A stat = 
(4.4 ± 0.6) x 10 3 s" 1 as quoted by Stocki et al [22]. 

We stress our goal was not to identify one particular model as preferable to another 
model. Rather, by different choices of valence nucleon spaces and single particle energies, 
we intended to expose the model dependences of the capture rate and the recoil alignment. In 
particular, by varying the ls-0<i single particle energies we changed the (0p)" 4 (ls0d) 2 wave- 
function admixture, the most obvious uncertainty in the nuclear structure of the 2 + doublet. 4 
This permitted a better understanding of how the uncertainties in the (0p)~ 4 (ls0d) 2 admix- 
ture are mapped into the uncertainties in the capture rate A stat and the recoil alignment 

j^stat 

A comparison of the excitation energies of the low-lying (J n , T) = (2 + , 1) states from the 
various models with the experimental data is given in Fig. 2. For CKPOT the single 2 + state 
is very nearly pure (lp 3 / 2 ) 7 (lpi/2) 3 - For MK3CW1 the lower-lying 2 + state is dominantly 
(lp 3 / 2 ) 7 (lpi/2) 3 and the higher-lying 2 + state is dominantly (0p)~ 4 (ls0d) 2 . Neither CKPOT 
nor MK3CW1 are capable of reproducing either the small energy gap of the 2 + doublet or 
the well-established mixing of the two 2 + states. However, for MK3CW2 and MK3CW3, in 
which we adjusted by hand the splitting between the Op shell and the IsOd shell, the two 
lowest-lying 2 + states are increasingly mixed, and the energy splitting of the 2 + doublet 
is thus adequately reproduced. Note for MK3CW2 the (0p)~ 4 (ls0d) 2 admixture is ~50%, 
this admixture dictating the choice of the single particle energies, while for MK3CW2 the 
(0p)~ 4 (ls0d) 2 admixture is ~85%, this admixture reproducing the capture rate A stat = 
4.4 x 10 3 s _1 . Clearly, the energy splitting of the 2 + doublet, the (0p) _4 (ls0d) 2 content of 
the 2 + states, and the relative positions of the Op, IsOd shells, are very strongly correlated. 

In order to compute the capture rate and recoil alignment for 14 N(1 + , 0)— » 14 C(2 + , 7012) 



4 The dependence of the (0p) _4 (ls0d) 2 wavefunction admixture on the ls-0d single particle en- 
ergies is a result of the proximity of the lowest-lying (Op) -2 and (0p) _4 (ls0d) 2 configurations 
with spin-parity 2 + . This circumstance implies a large mixing between (Op) -2 configurations and 
(0p) _4 (ls0d) 2 configurations in the 2 + doublet. 
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we employed the formalism of Walecka [55] For completeness we give in Appendix B the 
details of our model calculations of the recoil alignment. We took the weak nuclear current 
as a sum of A one-body nucleonic currents, i.e. ignoring the effects of exchange currents, 
isobar excitations, etc. The required nuclear matrix elements were computed with one-body 
transition densities obtained from the OXBASH shell model code [56] and with harmonic 
oscillator radial wavefunctions and an oscillator parameter b = 1.7 fm. We set the weak 
vector coupling g v — 1 and the weak magnetic coupling qm = 3.706, and assumed the 
induced scalar coupling g s and induced tensor couplings g t were zero. In order to scale 
the weak vector and magnetic couplings to finite q 2 we assumed a dipole form factor with 
A 2 = 0.73 GeV 2 . The three-momentum transfer q = 97 MeV/c was computed via energy- 
momentum conservation 

q + q 2 /2M t = - AE - e b (3) 

where M t is the target mass, AE is the 14 N— 14 C* nuclear binding energy difference, and 
€b is the muon binding energy. The muon wavefunction was assumed to be constant over 
the nucleus, with a reduction factor R = 0.84 accounting for the difference between a point 
nucleus and a finite nucleus (for details see Ref. [55]). 

Lastly, to assist our later discussions of model dependences, we make a few comments on 
the model calculations. We remind the reader a 1 + — > 2 + transition involves three multipoles 
jtt = 1+j 2 + ; 3+ ; although for 14 N(1+, 0)-> 14 C(2+, 7012) the contribution from J w = 3+ terms 
are very small. The J n = 1 + , 2 + multipoles involve four independent amplitudes (C\ — M\, 
j-eib _ j-mag ^ ^ _ j-ei _ j-magb^ ^adi themselves are products of basic multipole 
operators and weak coupling constants. For J 71 " = 1 + , 2 + the multipole operators are the 
allowed Gamow- Teller operator Mio • <x, second-forbidden Gamow- Teller operators M12 • u 
and M 2 2-o", and momentum-dependent contributions Mn-V, M 2 i-V, M 2 3- V and M^- V. 
The capture rate, recoil alignment and other observables in 14 N(1 + , 0)— > 14 C(2 + , 7012) are 
all combinations of the nuclear matrix elements of these basic multiple operators and the 
weak coupling constants. For further details see Walecka [55], Donnelly and Haxton [57] 
and Appendix B. 

C. Capture rate results 

The yU _ -|- 14 N(l + , 0) — > z/^+ 14 C(l + , 7012) capture rate for a statistical mixture of the 
hyperfine states, has been measured by Babaev et al. [39], Thompson et al. [40], Belotti 
et al. [41], G iff on et al. [42] and Stocki et al [22]. Unfortunately, the experimental results 
have significant scatter, with rates as small as (2.2 ± 0.9) x 10 3 s _1 [22] and as large as 
(10±3) x 10 3 s' 1 [39]. Stocki et al have quoted a world average of A stat = (4.4±0.6) x 10 3 s" 1 
[22]. 

Our results from the model calculations of the muon capture rate and its hyperfine depen- 
dence in the 14 N(1+, 0)^ 14 C(2 + , 7012) transition are listed in Table II. Also given in Table 
II are the (0p)~ 4 (ls0d) 2 admixture in the 14 C final state and the dominant 0p 3 / 2 — > 0pi/ 2 one- 
body transition density for the 14 N^ 14 C transition. Note that for 14 N(1+, 0) -> 14 C(2 + , 7012) 
the F = 3/2 capture rate and statistical capture rate are overwhelmingly determined by the 
product of the Gamow- Teller matrix element Mio ■ er and the weak axial coupling g a . More- 
over, as shown in Fig. 3, the capture rate is a strong function of the (0p)~ 4 (ls0d) 2 content of 
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the 14 C nuclear wavefunction, the rate reflecting the decreasing wavefunction overlap with 
increasing (0p)~ 4 (ls0d) 2 admixture. 

As found in earlier investigations the simple Op-shell grossly over-estimates the 14 N(1 + , 0) 
— > 14 C(1 + ,7012) capture rate. Clearly, increasing the 2hu admixture by simply changing 
the Op, ls-Od shell splitting is capable of eliminating the discrepancy between experiment 
and theory - although this procedure is rather ad hoc. 



D. Recoil alignment results 

Our results for the model calculations of the recoil alignment in the 14 N(l + ,0) — > 
14 C(2 + ,7012) transition are plotted in Fig. 4. Note that Fig. 4 shows the recoil align- 
ment for a statistical mixture of the hyperfine states, as determined by experiment and 
discussed in Sec. III. The alignment Af at of the statistical mixture was obtained from the 
two alignments Af of the hyperfine states via 

Astat _ f+A£ + f-A L R ,.s 
- U + f-R (4) 

where / + , /_ are the statistical populations of the hyperfine states and R = A_/A + is the 
hyperfine dependence of the capture rate. Therefore, in principle, the alignment A s ^ at is a 
function of the hyperfine alignments Af and the capture rates A ± , although, in practice, as 
R is small one finds the alignment is approximately A s ^ at ~ A\. 

A useful comparison for the model calculations of the recoil alignment is the Fujii- 
Primakoff approximation [58]. In the Fujii-Primakoff approximation the M 10 ■ u matrix 
element is retained while the forbidden matrix elements are omitted. The approximation 
yields 

A\ = +v / l4(-X 2 -X + 2)/(+4 + 2X + 5) (5) 




where the quantity X = (Ga — Gp)/Ga and the effective coupling constants are Ga = 
-(g a + (q/2M)(g v + g m )) and G P = -(q/2M)(g p -g a + g v + g m ). With g v = 1.00, g m = 3.706 
and q ~ 100 MeV/c one obtains a value of A\ ~ +0.12 for g p /g a — and A\ ~ +0.35 for 
g p / ' g a ~ 7, thus indicating the importance of the induced coupling on the recoil alignment. 

For g p /g a = 6.7 the model results of Fig. 4 give recoil alignments that range from Af at = 
+0.45 (MK3CW3) to Af at = +0.53 (CKPOT). The relatively small model dependence for 
the alignment may be contrasted with the relatively large model dependence of the capture 
rate, it decreasing by roughly a factor of five from CKPOT to MK3CW3 (see Table II). 
To better understand this striking difference in model sensitivities between A s ^ at and A stat 
we remind the reader the alignment is governed by ratios of nuclear matrix elements and 
weak coupling constants. Therefore, as both the numerator and denominator in A s ^ at are 
dominated by Mio • a, the recoil alignment A s ^ at is rather weakly model dependent although 
the matrix element M 10 • a itself is rather strongly model dependent. The residual model 
dependence of the recoil alignment A s ^ at may be ultimately traced to the small contributions 
of the ls-Od one-body transitions densities to the 14 N(l + ,0) — > 14 C(2 + ,7012) transition. 
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In Fig. 5 we show the relative contributions of the different matrix elements to the recoil 
alignment A s ^ at . The plot reveals a substantial contribution from the £ = 2 Gamow- Teller 
matrix elements (e.g. M 2 2 • cr) and a significant contribution from the axial charge matrix 
element (i.e. M^cr • V). The effects of other J = 2 matrix elements are small and the 
effects of the J = 3 matrix elements are negligible. Consequently, the dominant source of 
model dependences in computing A s ^ at is the knowledge of the ratio between the £ = 2 GT 
matrix elements and the Mio • u matrix element. This ratio is sensitive to the (0p)~ A (ls0d) 2 
wavefunction admixture through the ls-Od one-body transition densities. 

E. Coupling constant sensitivity 

Our alignment measurement, A s ^ at = 0.60±0.11, and alignment calculation, A s [ at = 0.45- 
0.53, are clearly in reasonable agreement. Indeed the comparison between model and data is 
sufficient to constrain the coupling's value to g p /g a > 5 (CKPOT) to g p /g a > 9 (MK3CW3). 
Our bound on g p /g a is consistent with the results from the 12 C recoil polarization experi- 
ments, which yielded g p /g a = 9.8 ± 1.8, but inconsistent with the results from the 28 Si recoil 
alignment experiments, which yielded g p /g a = l-Olj;^- Our result is also consistent with the 
prediction g p /g a = 6.7 of chiral symmetry arguments [2,3]. 

In considering the dependence of A s ^ at on g p its helpful to summarize the terms that 
contribute to A s ^ at (for definiteness we employ the MK3CW3 model). In the limit of (i) 
no forbidden nuclear matrix elements and (ii) no induced pseudoscalar coupling, the recoil 
alignment is Af; at ~ 0.09. The contribution of the coupling g p = 6.7g a increases the recoil 
alignment to A^ at ~ 0.25 and the contribution of £ = 2 GT terms increase the recoil 
alignment to A s ^ at ~ 0.40. Finally the remaining nuclear matrix elements, mostly the axial 
charge matrix element, are responsible for increasing the value to A s ^ at ~ 0.45. Clearly the 
agreement of model and data for A s ^ at lend support to significant contributions from both 
the induced pseudoscalar coupling and the £ = 2 forbidden corrections. 

V. SUMMARY 

In summary, we report a measurement of the longitudinal alignment Af at of the recoil 
nucleus in the 14 N(l + ,0) — > 14 C(2 + , 7012) transition. The alignment was determined from 
the Doppler lineshape of the 7012 keV gamma-ray from the 14 C(2 + , 7012) ^ 14 C(0 + ,0) 
transition, i.e. using the method of Grenacs et al. [23]. Our result of Af; at = 0.60 ±0.11 
indicates a large, positive alignment of the 14 C recoil nucleus originating from the V-A 
structure of the weak interaction. 

In addition, we report new calculations of the capture rate and the recoil alignment for 
the 14 N(1 + , 0) ^ 14 C(2 + , 7012) transition. The calculations were conducted in both a simple 
Qp-2 moc [ e i S p ace anc l a richer 0p~ 2 + 0p~ 4 (ls0d) 2 model space. As discussed by others, the 
capture rate A stat is dominated by the £ = Gamow- Teller matrix element and the weak 
axial coupling g a , and is highly sensitive to the (0p)~ 4 (ls0d) 2 admixture in the 12 C(2 + , 7012) 
wavefunction. By comparison, we found the recoil alignment A s ^ at has substantial contribu- 
tions from the induced pseudoscalar coupling and the second-forbidden Gamow- Teller matrix 
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elements. Interestingly, the recoil alignment is comparatively insensitive to the (0p)~ 4 (ls0d) 2 
components of the 12 C(2 + , 7012) wavefunction. 

Our measured alignment of A s ^ at = 0.60 ±0.11 and calculated alignments of A s ^ at = 
0.45-0.53 are in reasonable agreement. We view this reasonable agreement as evidence 
supporting the model estimates for the contributions of the induced pseudoscalar coupling 
and the second-forbidden corrections to the 14 C recoil alignment. Indeed, by comparing 
the model calculations and experimental data, we derived a bound on g p that ranges from 
9p/g a > 9 for the MK3CW3 calculation to g p /g a > 5 for the CKPOT calculation. This 
bound is consistent with the value g v j g a = 9.8 ±1.8 obtained from the determination of 
the 12 C(0 + ,0) ^ 12 B(l + ,0) transition recoil polarization, but inconsistent with the value 
g p /g a = l.Otil obtained from the determination of the 28 Si(0 + ,0) ^ 28 A1(1+, 2201) recoil 
alignment. Our result is consistent with g v j g a = 6.7 - the theoretical prediction derived 
from chiral symmetry arguments. 

Our experiment was limited by both the statistical accuracy of the Doppler lineshape 
and the unfortunate presence of the 35 Cl(n, 7) 6978 keV background line. The model uncer- 
tainties are apparently dominated by the uncertainties in the I = 2 Gamow- Teller matrix 
elements arising from the uncertainties in the (0p)~ A (ls0d) 2 wavefunction admixture. We 
suggest, with better statistics and reduced backgrounds, the 14 N(l + ,0) — > 14 C(2 + ,7012) 
transition would offer a competitive window on induced currents and forbidden contribu- 
tions in muon capture. 

We wish to thank Drs. Ermias Gete, David Measday, Belal Moftah and Michael Saliba 
for assistance in collecting the data and the TRIUMF technical staff for the operation of 
the TRIUMF cyclotron. We also acknowledge both the National Science Foundation (USA) 
and the Natural Sciences and Engineering Research Council (Canada) for their financial 
support. 
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APPENDIX A: RELATION BETWEEN THE 7-1/ CORRELATION AND THE 

LONGITUDINAL ALIGNMENT 



Herein we consider the sequence /i"+ 14 N(l + , 0) -> z/ M + 14 C(2+, 7012) -> 7 + 14 C(0 + ,0) 
which involves an allowed Gamow- Teller (1 + — > 2 + ) transition and pure electric quadrupole 
(2 + — > + ) decay. We denote the initial, intermediate and final state angular momenta by 
Ji, J2 and J3, the neutrino and photon directions by z> and k, and the 7-ray multipolarity 
by L w . For either unpolarized muons, P M = 0, or perpendicular geometry, P M • k = 0, the 
7-1/ directional correlation is given by [26-29] 

W{6) = l + J2 a s p s( cos9 ) ( A1 ) 

s 

where P s (cos9) are the Legendre polynomials, a s are the 7-1/ correlation coefficients, and 
cos 6* = u- k. In Eqn. Al the summation involves even Legendre polynomials, and is bounded 
by the smaller of 2J 2 or 2L. Thus for 14 N(1+, 0) ^ 14 C(2+, 7012) ^ 14 C(0+, 0) the permissible 
values are s = 2 and s = 4. 

The 7-z/ correlation coefficients a s are functions of both the \i capture process and 7-decay 
process. More specifically the coefficients a s are given by 



= A S B S (A2) 



where A s is dictated by the /i"+ 14 N(l+, g.s.) -> z/^+ 14 C(2+, 7012) transition, and B s is 
dictated by the 14 C(2+,7012) -> 7+ 14 C(0+,0) transition. The quantity A s is the rank-s 
orientation of the recoil nucleus J2 about the neutrino axis v. The quantity B s is governed 
by the spin-parity sequence J 2 — > J3 and the multipolarity L 7 " (a useful reference that 
contains tabulations of B s coefficients is Ciechanowicz and Oziewicz [29]). For 2 + — > + 

transitions one finds that B2 = — \J~^ and B4 = — ^/| . 

Our model calculations for /i~+ 14 N(l + , 0) -> z/ M + 14 C(2+, 7012) -> 7 + 14 C(0+,0) have 
shown the role of the rank-4 orientation in the 7-z/ correlation may be safely ignored. For 
the rank-2 orientation we have herein denoted A 2 as the recoil's longitudinal alignment 
Al = A 2 . The relationship between the 7-z/ directional correlation coefficient a 2 and the 

recoil longitudinal alignment A L is thus simply a 2 = —Jj^Al. 



APPENDIX B: MODEL CALCULATIONS OF THE LONGITUDINAL 

ALIGNMENT 



One approach in addressing the various physical observables in nuclear muon capture 
is the helicity representation. In the helicity representation one views capture as the two- 
body decay of a spin-F muonic atom into a left-handed muon neutrino and a spin- Jf recoil 
nucleus. The rank-s orientation Af of the recoil nucleus about the neutrino axis is 

Af = (-irv/27+I^T < J f \s0\J f \ > \7f | 2 / I^aT (Bl) 

A A 

where F is the muonic atom hyperfine state, A is the recoil nucleus helicity state, and T[ 
are the contributing helicity amplitudes [58,55,29]. 
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For /i"+ 14 N(l+,0) -> z/ M + 14 C(2+,7012) the initial hyperfine states are F = J t ± 1/2 = 
1/2,3/2. Thus according to Eqn. Bl, the various observables in F_ = 1/2 capture 
are determined by two helicity amplitudes, i.e. T^{ 2 , Tq 1 ^ 2 , and in F + = 3/2 capture 
are determined by four helicity amplitudes, i.e. T^ 2 , T^{ 2 , Tg^ 2 , T+{ 2 . These six ampli- 
tudes contain all the dependences on the weak couplings and the nuclear structure in the 
/L( -+ 14 N(1 + , 0) -> z/ M + 14 C(2+,7012) transition. 

More convenient for model calculations than the helicity amplitudes are the elec- 
troweak amplitudes denoted Cj — M.j and Tf — T™ a9 where A4j, £j, Tf and Tf ag rep- 
resent the charge, longitudinal, transverse electric and transverse magnetic operators. For 
/i"+ 14 N(l + , g.s.) — > z/ M + 14 C(2 + , 7012) there are six electroweak amplitudes corresponding 
to the three multipolarities J 71 " = 1 + , 2 + , 3 + where Cj—Mj and ^ji—^J 11 ^ are simple prod- 
ucts of basic multipole operators and weak coupling constants [55,57] (the specific relations 
between the electroweak amplitudes and the helicity amplitudes are given in Refs. [29]). For 
/i-+ 14 N(l+,0) -> z/ M + 14 C(2+,7012) the relevant formulas for Cj-Mj and Tf-T™ 9 with 
multipolarities J w = l + ,2 + are reproduced in Table III. They involve the allowed Gamow- 
Teller operator M 10 • a, the second forbidden Gamow- Teller operators M 12 ■ a and M 2 2 • o", 
and the momentum-dependent operators Mi cr • V, Mn ■ V, M21 • V and M 2 3 ■ V. The 
contributions from J w = 3 + terms to yU~+ 14 N(l + , 0) — > z/ M + 14 C(2 + , 7012) are very small. 

Finally, the required multi-particle weak matrix elements < J/||0 J ||Jj >, between an 
initial many-body state | Jj > and final many-body state \Jf >, were obtained from single- 
particle weak matrix elements < a'||0 J ||Q! >, between an initial single-particle state \a >= 
\n,j,£ > and final single-particle state \a' >= \n',j',£' >, via [57] 

< J f \\0 J \\Ji >= ^C(J,a,c/, J f , Ji) < a'||0 J ||a > (B2) 

a,a' 

where J denotes the relevant operator and C(J,a,a', Jf, Ji) denotes the so-called one- 
body transition densities. The one-body transition densities govern the contributions of the 
various single-particle matrix element < a'||0 J ||a > to the multi-particle matrix element 

<J f \\o J u>. 
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TABLES 



TABLE I. Summary of fits to the 7012 keV 7-ray lineshape from the 14 N(0, 1+) -► 14 C(7012, 2+) 
transition. Ac is the amplitude of the 6978 keV background peak, X D is the position of the 
7012 KeV Doppler peak, a is the instrumental resolution, and r s d is the slowing-down time. We 
also distinguish between the fits with a flat continuum background (denoted 1 par.) and a linear 
continuum background (denoted 2 par.). The resulting "best-fit" values for A s ^ at are listed in the 
last column. 



A G 


X D 


a G/D 


Tsd 


AE 


Bkgd. 




free 


free 


0.82 


0.011 


0.00747 


2 par. 


0.604 ±0.106 


free 


7012 


0.82 


0.011 


0.00747 


2 par. 


0.610 ±0.092 


free 


free 


1.08 


0.011 


0.00747 


2 par. 


0.617 ±0.115 


free 


free 


0.66 


0.011 


0.00747 


2 par. 


0.592 ± 0.095 


free 


free 


0.82 


0.000 


0.00747 


2 par. 


0.632 ±0.104 


free 


free 


0.82 


0.022 


0.00747 


2 par. 


0.578 ±0.105 


free 


free 


0.82 


0.011 


0.00747 


1 par. 


0.605 ±0.100 





free 


0.82 


0.011 


0.00747 


2 par. 


0.426 ± 0.078 


41 


free 


0.82 


0.011 


0.00747 


2 par. 


0.503 ± 0.083 
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TABLE II. Results of the model calculations for the capture rate and its hyperfine dependence 
for the 14 N(1 + , 0) — > 14 C(2 + , 7012) transition. The hyperfine capture rates are denoted A 1 * 1 and the 
statistical capture rate is denoted A stat . Also given for the different models are the 0p 3 / 2 — ► Opi/2 
one-body transition density (column five) and the 0p~ A ls0d 2 wavefunction admixture (column six). 



Model 


A+ 


A_ 


^stat 


A-/A+ 


0P3/2 QPl/2 


4h-2p 


CKPOT 


24.8 x 10 3 


1.0 x 10 3 


16.9 x 10 3 


0.040 


1.100 


0% 


MK3CW1 


23.8 x 10 3 


1.3 x 10 3 


16.3 x 10 3 


0.053 


1.005 


~20% 


MK3CW2 


18.3 x 10 3 


1.1 x 10 3 


12.6 x 10 3 


0.058 


0.808 


~50% 


MK3CW3 


6.3 x 10 3 


0.42 x 10 3 


4.3 x 10 3 


0.067 


0.362 


-85% 



TABLE III. Decomposition of the electroweak operators Cj — Aij and Tj — jj 1 ^ w ith multi- 
polarities J 71 " = l + ,2 + into the basic multipole operators and the weak coupling constants. The 
basic multipole operators Mj, Mjl ■ <r, Mjl ■ V and Mj a ■ V are discussed in the Appendix B. 





electroweak 


multipole operator 






operator 


decomposition 




1+ 


L\ - M\ 


i[(9a + ^(Oa ~ g p ))(\[l M w ■ a + y/lM 12 ■ tr) 


+ ±g a M 1 cr • V]r± 


1+ 


rpelb _ rpmag 


i[(ga - m(9v - g m ))(\flM 10 ■ a - \J\M 12 ■ a 


) + favMu • V]r± 


2+ 


L 2 - M 2 






2+ 


rpel _ rpinagh 


i[{-9a + £l(9v + g m ))M 2 2 ■ <r + ±g v {JlM 2 i ■ 


V-TfM^-V)]^ 
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FIGURES 
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FIG. 1. The Doppler lineshape of the 7012 keV gamma-ray from the fi + N(l + ,0) — > 
i/ M + 14 C(2+,7012) -► 7+ 14 C(0+,0) sequence. The rounded top of the 7012 keV lineshape is a 
reflection of the 14 C recoil alignment in the \i capture process. The small peak on the far lefthand 
side of the Doppler lineshape is the 6978 keV background line. The solid line is the benchmark fit 
to the Doppler lineshape (see text for details). 
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FIG. 2. Energy levels of the low-lying 2 + states from the experimental data and the model 
calculations. The experimental excitation energies of the first and second 2 + states are respec- 
tively 7012 and 8318 keV. The 2fiu admixture in the lowest-lying 2+ state of the three MK3CW 
calculations is increasing from left to right. 
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100 



1s0d admixture (%) 

FIG. 3. The statistical capture rate for ^"+ 14 N(1 + , g.s.) -> i/ M + 14 C(2+, 7012) as a function of 
the (0p)~ A (ls0d) 2 wavefunction admixture. The calculations were performed by varying the energy 
separation between the Op shell and the IsOd shell in the MK3CW model calculations. 
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4 8 12 16 20 

gp/ga 

FIG. 4. The recoil alignment of the 14 C recoil as a function of the coupling constant ratio g p /g a - 
The calculations were performed for CKPOT (diamonds), MK3CW1 (open circles), MK3CW2 
(filled circles), and MK3CW3 (stars). The horizontal lines indicate the central value and error bar 
for the experimental result A s ^ at = 0.60 ± 0.11. 
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4 8 12 16 20 

gp/ga 

FIG. 5. The recoil alignment of the 14 C recoil as a function of the coupling constant ratio g p /g a - 
The calculations were performed with MK3CW3 for M\Q-a only (diamonds), for M\q-g and Mj^-cr 
(open circles), for M\q ■ a, Mj2 ■ <r and M±a ■ V (filled circles), and with for J = 1,2,3 matrix 
elements (stars). The horizontal lines indicate the central value and error bar for the experimental 
result A s l at = 0.60 ± 0.11. 
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